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Abstract 

New types of symmetry for the Rayleigh equation are found. For small Atwood num- 
ber, an analytic solution is obtained for a smoothly varying density profile. It is shown 
that a transition layer with a finite width can undergo some kind of stratification. 
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Although Rayleigh had investigated the instability (known now as Rayleigh- Taylor (RT) 
instability) using an exponential profile of density most investigations have been done so 
using sharp interfaces [2-6]. In [0] it is carried out making an interesting analogy between the 
equations which describes the RT instability, and the Schrodinger equation. In that study, an 
"equivalent potential" is constructed and for well-known density profiles it is shown that the 
Schrodinger equation with the corresponding "equivalent potential" has the same eigenvalues 
as the Rayleigh equation. 

The physical quantities (e.g.; density, velocity) and their derivatives, generally speaking 
cannot suffer a jump discontinuity. Therefore, it seems more consistent to attempt to solve the 
problem for a transition layer of finite thickness and then take the limit when the thickness of 
this layer A tends to zero; i.e. consider the case of a density jump in the limit, A — ► 0. The 
solution for this transition layer gives us a more complete physical picture of the instability. 

In this paper we consider the RT instability of a transition layer of finite thickness, where 
the unperturbed density changes continuously from a constant value up to another one. For 
that, an analytic solution is found, in the limit of small Atwood number. This solution allows 
us to investigate the A — ► limit. We find that: 

1. For an arbitrary density profile the "equivalent potential" can be presented in compact 
form. This allowed us to predict a new type of symmetry (in addition to the well-known 
symmetry of the Rayleigh equation @). 

2. The depth of the "equivalent potential" well depends on the width of the transition layer 
and the depth not the energy of the potential well undergoes a quantization. For the finite 
width of the layer the spectrum of eigenvalues is infinite. This conclusion is consistent with the 
fact that with the increase of the quantum number grows the well depth. 

3. The eigenf unctions found here show the exfoliation of the transition layer. 

I. Using the linearized equation of motion and the continuity equation with the help of 
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incompressibility condition leads one to: 

where po(z) describes the equilibrium density profile, v(z, t) is the z component of the perturbed 
fluid velocity, k is the wave number of the perturbed quantities which are chosen proportional 
to e lky and g is the constant acceleration along negative z. Then one can find the solution of 
the above equation using separation of variables: 

v(z,t)=T(t)$(z). (2) 

Then Eq. (1) reduces into two equations as follows: 

T"(t) - 7 2 T(t) = 0, (3) 

+ ^t4^) + (~k 2 + ^PM)$( Z ) = 0, (4) 
PoO) T Po(z) 

where 7 2 is the constant of separation, primes denote derivatives with respect to corresponding 
arguments. The case 7 2 > corresponds to unstable modes. 
The substitution 

reduces Eq. (4) to the form of the Schrodinger equation as in [0] : 

i)"(z) + [-k 2 - V(z)]ip(z) = 0, (6) 

with the "equivalent potential" , 
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where 

r kg 

7 2 ' 

One can further simplify Eqs. (6) , (7) as follows: 

S"(z 



V(z) = ^ - Pr 2 , (9) 

1/)"(Z) , 2 S"'(z) , 2 . 

ip(z) S(z) 



where 

S(z) = e- krz ^p- . (11) 

This form then allows us to prove the following symmetry properties of the Rayleigh equation 
for arbitrary density profiles: 

a) The substitutions 

T^-r and p (z) - p (z)e- 4krz (12) 

cause Eq. (9), and consequently Eq. (6), to remain unchanged. This means that every unstable 
mode, T > 0, with the density profile po(z) can be compared with a stable, oscillating mode, 
T < (see Eqs. (3) , (8)), with the density profile po(z)e~ 4kTz and vice versa. For the modified 
profile, on the other hand, the expression for the velocity (see Eqs. (2) , (5)) will contain an 
additional factor e 2kTz . 

b) Performing the substitution po(z) by po(— z)e AkTz and then the transformation z — > —z, 
we find that the potential (9) and Eq. (6) remain unchanged. Consequently the spectrum of the 
eigenvalues for both profiles are the same, whereas eigenfunction of the latter can be obtained by 
sending z — > —z in the former. Furthermore, the function $(z) upon transformation becomes 
e~ 2kTz &(z). Note that these symmetry properties hold only for fixed Y n (n = 0, 1, 2, ...) 

II. Let us assume that the continuous equilibrium density tends to different constant limits, 
i.e. po(z) — > p± at z — > — oo and po(z) — > p2 at z — > +oo, where P2 > Pi , and p' (±oo) = 0. 
Then if somewhere in the region — oo < z < +oo the "equivalent potential" becomes negative, 
according to Eq. (6) the RT problem becomes analogous to a particle's motion in a potential 
well. Using a tanh function representation for density: 

2z 

p (z) =c(l+A-tanh— ), (13) 

where 

c = Pi+Pi A = ^ (14) 

2 P2 + Pi 

enable us to examine different limits, from an smooth to a sharp jump limit. We restrict 
ourselves to the case of small Atwood number , A < 1, where one can obtain an analytic 



solution. Neglecting the terms of order A 2 from Eq. (6) we obtain: 

i//'(x) + {-<5 2 + — -g-[<J ■ AT + A ■ tanhx]}^(x) = 0, (15) 
cosh x 

where 5 = kA/2 , x = 2z/A. We made an assumption that AT > 1, which will be confirmed 
by the results. For a broad transition layer, i.e. 5 ^> A, Eq. (15) can be reduced to a form 
which is well-known in quantum mechanics M (see p. 73): 

x . Ar 

r{x) + H 2 + —^Mx) = (16) 
cosh x 

Solution to this equation finite for z — > ±oo, can be expressed in terms of hypergeometric 
function F: 

ip = cosh -6 x F(6- s,6 + s + l,5 + l,(l - tanh x) /2), (17) 



the parameters of which (5 , s = |[— 1 + yl + 45 AT] ) have to satisfy the relation: 

«5- s = -n, (18) 

where n — 0, 1, 2, ... . From Eq. (18) we find the spectrum of growth rates for RT instability: 

7 2 5 

— = A (19) 

kg (6 + n)(5 + n + l) 

Consequently for a transition layer of finite width, 5 ~ 1, there is an infinite spectrum of 
eigenvalues. Number of the extremum and zero points of the eigenf unctions (17) is equal to 
n + 1 and n , respectively. For n > the spatial dependence of the eigenfunction has oscillatory 
character. It means that due to the instability, the monotonous distribution of the density in 
the transition layer can be destroyed, and a stratified and multi-stream structure can be formed. 
The velocity changes the direction as one goes from one new-formed sublayer to another. The 
Figure shows the spatial distribution of the velocity (determined from the function <&(z); see 
Eq. (2)) and the total density p(z) = po(z) + p(z) in such a stratified layer. The density 
perturbation p(z) is obtained from the continuity equation: 



As we are concerned in the spatial dependence only, by plotting the figure we treated the 
definite integral as constant (for simplicity we choose it as equal to one). If the width of the 
transition layer is small, say 

5 = kA/2 <A<1, (21) 

we can consider such a layer as a jump in density. However, in this case the analytic solution 
of Eq. (15) fails. It is necessary to note that although the spectrum represented by Eq. (19) is 
obtained for a broad transition layer, 5 > A, (when we neglect the second term in the square 
brackets in Eq. (15)) Eq. (19) gives the correct expression for the main unstable mode in the 
case of density jump, when 5 < In fact as 5 — > from Eq. (19) it follows that the most 
unstable mode corresponds to the quantum number n — 0, and for the growth rate one obtains 
the well-known expression JTJ. 0: 

7 = \fkgA- (22) 

Such a general meaning of Eq. (19) can be explained by the fact that, for 5 C 1 , A <C 1 
the "equivalent potential" in Eq. (15) can be considered as a perturbation. The quantum 
mechanical problem of particle motion in a one-dimensional well, whose depth is small, is 
solved in [|J (see p. 162). There an equation for the determination of the eigenvalues of the 
Schrodinger equation is obtained. That equation in our notation becomes: 

/+oo 
dxV(x). (23) 
-oo 

It is obvious that the second term of Eq. (15) doesn't give any contribution in this integral 
(this term is odd). The first term gives the correct expression for the growth rate consistent 
with Eq. (22). That is why for the small "equivalent potential", (i5 < 1,A C 1) the second 
term in square brackets in Eq. (15) does not play any role in the determination of the growth 
rate and the expression (19) is valid for a shallow well also. 
We can summarize our results as follows: 

1. The new symmetry properties found in this paper show new density profiles for which 
the problem of instability can be reduced to the one and the same equation 0, [7|, [!|. 
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2. According to the symmetry theorem known up to now, the two profiles po(z) and l/po(— z) 
have the same set of eigenvalues but not eigenfunctions [0, ||, |9[]. The symmetry properties found 
here maintain the spectrum of eigenvalues and allow one to determine a relationship between 
the eigenfunctions (see Eq. (5)). 

3. The Eq. (10) shows that the inverse problem (obtaining the density profile for known 
perturbed quantity) is described by the same type of equation as the original problem. 

4. From Eq. (17) it follows that for every eigenmode, there exists a unique transition layer 
of stratified structure. 

5. It should be interesting to extend this analysis when separation of variables is not used 
and thus the spatial and temporal parts not treated independently. 

These results are true not only for the chosen density profile from Eq. (13) but also for any 
density profile whose corresponding "equivalent potential" has the form of a well. 

We are indebted to prof. N.L.Tsintsadze for encouragement and discussions. We thank Dr. 
F. Kazeminezhad for help in preparing of manuscript. 
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